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Abstract. These lectures describe some aspects of the physics of massive neutrinos. 
After a brief introduction of neutrinos in the Standard Model, I discuss possible pat- 
2 ' terns for their masses. In particular, I show how the presence of a large Majorana mass 

term for the right-handed neutrinos can engender tiny neutrino masses for the observed 
iy~j ' neutrinos. If neutrinos have mass, different flavors of neutrinos can oscillate into one 

04 , another. To analyze this phenomena, I develop the relevant formalism for neutrino 

oscillations, both in vacuum and in matter. After reviewing the existing (negative) 
evidence for neutrino masses coming from direct searches, I discuss evidence for, and 
hints of, neutrino oscillations in the atmosphere, the sun, and at accelerators. Some of 
^^ . the theoretical implications of these results are emphasized. I close these lectures by 

lO ' briefly outlining future experiments which will shed further light on atmospheric, accel- 

>0 , erator and solar neutrino oscillations. A pedagogical discussion of Dirac and Majorana 

^r^ ' masses is contained in an appendix. 

-p,; I NEUTRINOS IN THE STANDARD MODEL 

qj I Neutrinos play a special role in the SU{2) x f/(l) electroweak theory. While the 

^ ■ left-handed neutrinos are part of SU{2) doublets 

><: L,= (^'^»^ , £, = {e,/i,r}, (1) 

the right-handed neutrinos are SU{2) singlets. Since the electromagnetic charge 
and the ^(1) hypercharge differ by the value of the third component of weak isospin 

Q = T, + Y, (2) 

one sees that the right-handed neutrinos (t'^Jn carry no SU{2) x f/(l) quantum 
numbers. The above has two important experimental implications: 

i) The neutrinos seen experimentally are those produced by the weak interactions. 
These neutrinos are purely left-handed. 

ii) Because one cannot infer the existence of right-handed neutrinos from weak 
processes, the presence of (t'^Jn can only be seen indirectly, most likely through 
the existence of neutrino masses. However, it should be noted that neutrino 
masses do not necessarily imply the existence of right-handed neutrinos. 



Before discussing issues connected with neutrino masses, it is useful to summarize 
what we know about the left-handed neutrinos from their weak interactions. In the 
electroweak theory these neutrinos couple to the Z-boson in a universal fashion 

I^Zvv = ^ ^ ■ n Z^[J ]y , (3) 

where 

[J^,% = E(^^Jl7m(^4)l . (4) 

i 

Precision studies of the Z-line shape allow the determination of the number of 
neutrino species i, since as this number increases so does the Z total width. Each 
neutrino type (provided its mass m^^ <^ ^) contributes the same amount to the 
Z-width [1] 
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T{Z -^ u,M,^ = ,,1 ' p . (5) 



Here Gp is the Fermi constant as determined in /i-decay [2] and p is related to the 
axial coupling of the charged leptons to the Z-boson: p = (25f^)^. Using [1] 

^^ = -0.50102 ± 0.00030 , (6) 

which is the average of the results obtained by the four LEP collaborations and 
SLD, one has numerically 

r^ = T{Z -^ ueMeJ = (167.06 ± 0.22) MeV . (7) 

Using the above, the number of different neutrino species, N,y, can then be derived 
from the precision measurements of the Z total width and of its partial width into 
hadrons and leptons using the (obvious) equation 

Ttot = Thad + SFiept + N^T^, . (8) 

From a fit of the Z-line shape for e^e^ -^ p^ p^ and e^e^ -^ hadrons at LEP 
one can extract very accurate values for Ftot, Tiept and Fhad: 

Ftot = (2.4939 ± 0.0024) GeV 
Tiept = (83.90 ±0.1) MeV 
Thad = (1.7423 ± 0.0023) GeV . (9) 

Whence it follows that the, so-called, invisible width is 

Finv = N^V^ = (499.9 ± 3.4) MeV (10) 

and thus one deduces for N^, — the number of neutrino species: 



iV^ = 2.992 ± 0.020 . (11) 

This result strongly supports the notion, expressed in Eq. (1), that there are only 
3 generations of leptons. 

It turns out that one can get a more accurate value for Ni, by using other infor- 
mation derivable from the Z-line shape. The cross-section for e^e~ — > hadrons can 
be expressed in terms of three factors: [3] a peak cross-section 
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a Breit-Wigner factor 
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and a computable initial state bremsstrahlung correction (1 — 5qed(s))j with 

^had = ^oBW(s)(l-(5QED(s)) . (14) 

The value of o"o extracted from an analysis of the cross-section for e^e^ -^ hadrons 
at LEP [1] 

ao = (41.491 ± 0.058) nb (15) 

can be combined with the LEP results for the ratio of hadronic to leptonic partial 
widths of the Z 

R^ = -J^ = 20.765 ± 0.026 (16) 

1 lept 

to deduce, with a little bit of theoretical input, a value for N^. This value, as we 
will see, is slightly more accurate than that given in Eq. (11). 
Using Eqs. (8), (10), (12), and (16), one can write 



Tinv _ Tlcpt j Ttot Thad ol _ "'^lept j /I2vr-R 



^^^Mm.^Mcpt_^to^_M^_3lMept t^!l_^ - i?, - 3 . (17) 



In the Standard Model, the ratio Ficpt/Pj, is very accurately known: 



lept 



= 1.991 ±0.001 . 

SM 



Using this value in Eq. (17), along with the experimentally determined Z mass 
Mz = (91.1867 ±0.0021) GeV and the values of Uo and Ri measured at LEP, gives 

A^^ = 2.994 ± 0.011 ; 1;,^ = (500.1 ± 1.9) MeV . (19) 



These values are consistent with those in Eqs. (10) and (11), but are about a factor 
of two more accurate. 

There is an analogous equation to Eq. (3) describing the coupling of the W^ bo- 
son to the leptonic charged currents. Again only left-handed neutrinos are involved. 
One has 

Cwi., = ^ ' , {Wpi^' + W^Jif} , (20) 

V 2 sm Ow 

where 

J^^* = (^ST=E^"^L7.^4L. (21) 

i 

The states that appear in Eq. (21) in general are not mass eigenstates, since mass 
generation can mix leptons of the same charge among each other. Nevertheless, 
one can always diagonalize the charged lepton mass matrix by a by-unitary trans- 
formation of the left- and right-handed charged lepton fields: 

4 = U% ; in = V% . (22) 

After this transformation, the charged currents in Eq. (21) read 

f;^' = {fXy = E ^Ll,{U%ue^^ = E ^Ll,i^e.L , (23) 



y 



where 



^>£L = {UyueL . (24) 



Note that because U is unitary, the neutral current [J^ ],/ is the same whether 
it is expressed in terms of z/^l or z/^l- Conventionally, the states z/^-l are called 
weak interaction eigenstates, since they are produced in the decay of a W~^ 
boson in association with a physically charged lepton ii of definite mass. These 
states, of course, are also pair produced by the Z-boson. For ease of notation, in 
what follows I will drop the tilde on both £jl and z/^.^ with the understanding that 
the states now called z/^.l are those produced by the weak interactions — they are 
the weak interaction eigenstates. Similarly, the charged leptons ii are the states 
associated with a diagonal mass matrix 



Mf = \ m 



rUf. \ 



M 
m^ 



(25) 
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II PATTERNS OF NEUTRINO MASSES 

With these preUminaries underway, I want now to examine in a bit of detail 
the possible patterns of neutrino masses. To do so, it is useful to first review how 
fermion masses originate in field theory. The mass term with which everybody is 
acquainted with is one involving a fermion ip and its conjugate i/j = ip'^'-f^: 

^mass = -mi)i) = -m(V^L^R + ^r^l) , (26) 

with ipi^, i/jji being the usual projections 

V-L = ^(1 - 75)V' ; V-R = 1(1 + 75)V' • (27) 

This term, obviously, conserves fermion number 

ip -^ e^'^ilj ; i/j -^ e-'^i/j (28) 

and gives equal mass for particles and antiparticles 

m^ = m^ = m . (29) 

For particles carrying any f/(l) quantum number, like electromagnetic charge, it 
is clear that £mass is the only possible mass term, since to preserve these U{\) 
quantum numbers one needs always to have particle-antiparticle interactions. 

Neutrinos, however, provide an interesting exception. Because neutrinos do not 
have electromagnetic charge, it is possible to contemplate other types of mass terms 
for them besides the particle-antiparticle term given in Eq. (26). These other 
neutrino mass terms, contain two neutrino (or two antineutrino) fields. Hence they 
violate fermion number (and in some cases SU{2) x f/(l)), but otherwise are allowed 
by Lorentz invariance. 

As we discuss in more detail in Appendix A, one can write three different types 
of mass terms for neutrinos: 

- \ {^ICmTV^ + V^Cm}T:u^^\ . (30) 

Here the mass matrices vriD^msiVaT are Lorentz scalars. However, their presence 
is only possible as a result of different symmetry breakdowns. Specifically, vhd 
conserves fermion number, but violate SU{2) x U{1) since it does not transform 
as an SU{2) doublet. This fermion number conserving mass is often called a Dirac 
mass. Thus, in a happy confiuence of notation, m£, can stand both for a Dirac mass 
and a doublet mass. Both rris and m'r violate fermion number by two units and 
are known as Majorana masses. Because ms couples z/r with itself, clearly it is an 



SU{2) X f/(l) invariant. This is not the case for mrp, which violates SU{2) x f/(l) 
because it does not transform as an SU{2) triplet. 

The matrix C which enters in the Majorana mass terms in Eq. (30) is there 
to preserve Lorentz invariance. Appendix A contains a detailed discussion of this 
point, along with a pedagogical review of how one constructs 4-spinors starting from 
2-dimensional Weyl spinors. I note here only that C is not to be confused with the 
matrix C connected with how Dirac fields transform under charge conjugation. [4]^ 
Under the charge conjugation operator U{C) a Dirac field ip is transformed into its 
Hermitian conjugate ip"^: 

U{C)^U{C)-^ = C^\x) . (31) 

The matrix C is necessary to insure the invariance of the Dirac equation under 
charge conjugation and obeys the restriction 

Ci;C-' = -7, . (32) 

In general, C depends on the 7-matrix basis used. In the Majorana basis, where the 
7- matrices are purely imaginary, then C = 1. At any rate, the matrix C appearing 
in Eq. (30) is related to C by [5] 

C = CY^ . (33) 

It is easy to check that instead of (32) C obeys 

Clj:C~' = -7. . (34) 

The reason the matrix C appears in Eq. (30) is because it relates the, so-called, 
charge conjugate field ip'^ to ip rather than to ip'' . In view of the way the charge 
conjugation operator acts on the fermion field ip (see Eq. (31)), it is natural to 
define the charge conjugate field ip'^ as 

^^(x) = C^\x) . (35) 

Now ip = ip'^'j^, so one also has that 

^-(x) = C^^^^^ix) = C^^{x) . (36) 

So C, indeed, serves to relate ip to ip'^. 



^' Unfortunately, the distinction between C and C is often blurred in the literature. 



Ill THE SEE-SAW MECHANISM 

Eq. (30) displays the most general neutrino mass term, involving three distinct 
mass matrices mrt-, irts and itit- If there are only three flavors of neutrinos these 
are 3x3 matrices. For the moment this is what we shall assume, but we shall 
return to this point later on. 

One can write Eq. (30) in a more symmetrical way by replacing the transposed 
fields in this equation by the charge conjugate field. Recall that [cf. Eq. (35)] 

^^ = C^^ ; ^ = ^^C . (37) 

Hence, for example, one can write^ 

Wi^L = -i^lw^ = i^lCCw^ = Ulu^ . (38) 

Thus 

z^z^L = T^pK^h +'i4.^^ ■ (39) 

Using these equations, C^^^^s '^^^ be written in the following compact way: 

1 



r^ -_ 

mass r) 



f— \ ( rriT m^\( ^l 

^'- ^^) [mo ms [u^^ 



h.c. (40) 



For 3 generations of neutrinos, the six mass eigenstates rrii are the eigenvalues of 
the 6 X 6 matrix 

M=f ^^ <] . (41) 

yniD ms ) 

Because M is not necessarily Hermitian, its diagonalization necessitates a bi-unitary 
transformation 

UlMUi^ = Mdiag , (42) 

where J/l and f/a are 6x6 unitary matrices. This diagonalization is accomplished 
by a basis change on the original neutrino fields 

to a new set of fields ?7l and ?7r defined by the equations: 

^L = Ui^Vh ; ^R = f^R^R • (44) 



The minus sign in the second term below comes from Fermi statistics. 



It is useful to consider the simple, but physically interesting case, [6] of just one 
family of neutrinos. Further, let us imagine rriT = and ms ^ mo- The 2x2 
matrix M in this case reads simply 

M=f ^^"l . (45) 

This matrix has two eigenvalues, given approximately by ms and —mj^/ms- That 
is, in this case the spectrum splits into a very heavy neutrino of (approximate) mass 
m,s and a very light neutrino of (approximate) mass m^j~,/ms.^ This, so called, see- 
saw mechanism is very suggestive. It is natural to expect that m^ should be of 
the order of the charged lepton mass, corresponding to the neutrino in question: 
niD ~ mil. Then the spectrum of leptons has a natural hierarchy: 

("^/.)iight ~ m^ f — ^ j < m£ < (m^)hcavy ~ ms . (46) 

So, if there is a large mass scale associated with the right-handed neutrinos (the 
mass scale m,s, which is not constrained by the scale of SU{2) x U{1) breaking, 
since it is an SU{2) x f/(l) singlet) one readily understands why neutrino masses 
could be so much lighter than the corresponding charged lepton masses. 

The matrix M in the simple 2x2 example above is diagonalized (approximately) 
by the orthogonal matrix 

U=( \ rnn/rns] 

\ -mo/ms 1 / 

The two neutrino mass eigenstates are then 

"^"--i m ).^[ mo/ms 1 ' ][ ul ] (^^ 





I note that the mass eigenstates rji and 772 are Majorana (self-conjugate) states 

m = ^iL + ^iR = («^L + i^l) ^('^R + «^r) = Vi (50) 

ms 

V2 = ^2L + ^2R = (Z^R + I^r) H -il^L + l^l) = Vl ■ (51) 

ms 

The z/L state which enters in the weak interactions, for all practical purposes is, 
essentially rjn^. That is, it is the state associated with the light neutrino eigenstate 

[rrii ^ mjj/ms)'. 



^' For fermion fields, the sign of the mass term is irrelevant since it can be changed by a chiral 
transformation tpji -^ cxp [z^] ^r; i/)l — *■ exp [—if] V'L which leaves the rest of the Lagrangian 
invariant. 
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i^L = rjVL H ri2L ■ (52) 

ms 

The right-handed neutrino z/r, on the other hand, is essentially the heavy neutrinos 
eigenstate rj2R {jn2 — ms): 

^R = ri2K ^IR • (53) 

ms 

This simple example can be easily generalized to the 3x3 case of interest. Again, 
if the matrix itlt is negligible (i.e. if its eigenvalues are negligibly small), then the 
neutrino mass matrix M takes the approximate form 

M=f ^S") . (54) 

\mD ms ) 

Provided the eigenvalues of ms are large compared to those of ma, then again the 
spectrum separates into a light and heavy neutrino sector. The light neutrinos have 
a 3 X 3 mass matrix 

(Mi,)iight = m'^mg^mo , (55) 

while the heavy neutrino mass eigenstates are the eigenstates of the 3x3 matrix 

(M^)hcavy = ms . (56) 

The see-saw mechanism, in my view, is the only natural way to understand eV 
neutrino masses. Let me expand a bit on this point. Simce ms is an SU{2) x f/(l) 
invariant parameter, there are no constraints on it. On the other hand, as we 
discussed earlier, both m,£, and m-r can only originate after SU{2) x f/(l) breaking. 

The Yukawa interaction, of z/r with a left-handed doublet L = „ ] via a Higgs 
doublet $ = I ^! J 

-^Yukawa = -Fz^r^L + h.C. , (57) 

leads to a Dirac mass 

mo = r(0°) . (58) 

Since {(p^) is fixed by the scale of the SU{2) x f/(l) breakdown: 

(0°) = —^ ~ 180 GeV , (59) 



to get m,D to have a value in the eV range requires that F ~ 10^ 
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The situation is not much less artificial in the case of rriT- In this case, to get a 
non-zero value for ttit it is necessary to introduce a Higgs triplet field A. This field 
can couple to L L so that if, indeed, A gets a VEV one can generate a triplet 
mass rriT- In detail, the triplet coupling involving A has the form 

Aripiet = —{TtL^Ct ■ AL} + h.c. , (60) 

where T^ is an unknown coupling constant. When the neutral component of A, 
A°, gets a vacuum expectation, then triplet generates a mass term for z/l: 

Cass = -^r^(A°){z/?C'z/L} + h.c. (61) 

and rriT = Tt{A°). The only real constraint on (A°) comes from precision mea- 
surements of the p parameter, typifying the NC to CC ratio. Experimentally [1] 
one finds 

Pcxp = 1.00412 ± 0.00124 . (62) 

The presence of the triplet Higgs interaction modifies the p parameter from unity 
at the tree level and one has: [7] 

p = 1 - 2 [ -^— I J + rad. corr. (63) 

Using the error on p in Eq. (62) as an estimate of the size of (A°) implies that 
(A") < 4 GeV. So, also in this case, if (A°) is near this limit to get neutrino masses 
in the eV range one needs a Yukawa strength of order Ft ~ 10~^. If (A°) << {(p^) 
then Ft can be larger, but one is left to explain the reason for the doublet-triplet 
VEV hierarchy. 

Elementary Higgs triplets do not emerge very naturally in models. ^ However, 
one can always get an effective triplet out of two Higgs doublets: 

A ~ <I)^Cf$ (64) 

where C is an appropriate charge conjugation matrix. Effective L- violating interac- 
tions involving pairs of doublet Higgs fields arise quite naturally in Grand Unified 
Theories [9] as dimension 5 terms: 

>Cfr = ^(L^CrL) . {^^Cm + h.c. (65) 

In the above, A is a scale associated with the GUT breakdown scale and g is a 
coupling constant. Clearly the above interaction gives 



'^' An exception is provided by left-right symmetric models where triplets have often been con- 
sidered to give the requisite symmetry breaking. [8] 



rriT = ^^ . (66) 

Since (0°) ~ 10^ GeV, with g ~ 0(1), one gets itit ~ 10"^ eV for scales A ~ 10^^ 
GeV, which are typical of GUTs. Note that the above formula for rriT is quite 
similar in spirit to the see-saw expression for light neutrinos 

r,^ \see-saw ^^D \^ ) ( (-^\ 

since mci ~ (0°)- In either case, new physics at a large scale (either a large z/r 
mass scale ms or the GUT scale A) produces a light neutrino. It is clearly more 
appealing physically to have light neutrinos be the result of new physics at high 
scales, rather than simply as a result of some Yukawa coupling being unnaturally 
small. 



IV NEUTRINO OSCILLATIONS IN VACUUM 

If neutrinos have mass then, in general, the neutrinos produced by the weak 
interactions (weak interaction eigenstates) are not states of definite mass (mass 
eigenstates). In the basis where the charged lepton mass matrix is diagonal [c.f. 
Eq. (25)], it follows from Eq. (23) that the neutrino weak interaction eigenstates 
are fixed by the corresponding lepton produced in the associated weak process. 
That is, the piece of the weak current J^T involving the charged lepton £j = 
{e, /i, r} will always involve the corresponding neutrino weak interaction eigenstate 
Vi^ = {z/e, z/^, z/r}- These, left-handed, neutrino weak interaction eigenstates are 
superpositions of neutrino mass eigenstates z/j: 

'^'^j = J2 Uiji^i ■ (68) 

i 

The matrix f/^ j, in general, is a 3 x 6 matrix. However, if the see-saw mechanism 
is operative, one expects that the contributions of superheavy neutrinos in Eq. 
(68) should be negligible. Then, to a very good approximation, the matrix U^^i is 
a 3 X 3 unitary matrix 

U,^,Ul, = 5e^e, • (69) 

For the moment I will restrict myself to the case when Eq. (69) holds. Further- 
more, I will discuss the phenomenology of neutrino oscillations in the simple case 
of just 2 flavors of neutrinos, since this is how most of the data is usually presented. 
However, the formalism which we will develop can be generalized straightforwardly 
to three families of light neutrinos. [10] For deflnitiveness, let us consider then 
just z/g and z/^ weak interaction eigenstates. In this case, Eq. (68) reads, using a 
convenient quantum mechanical notation. 



|z/e) = cos^|z/i) + sin^|z/2) 
\u^) = — sin^|z/i) + cos^|z/2) • 



(70) 



The mass eigenstates |z/j) have a time evolution which just follows from the 
Schrodinger equation: 



W^it)) 



-iE^tl 



^^m 



E, 



p'^ + mf 



(71) 



Because mi ^ 1712, it is easy to see that the weak interaction eigenstate z/g 
produced at t = evolves in time into a superposition of z/g and u^ states. Taking 
by definition |z/j) = |z/j(0)), it follows that 



|z/e(t)) = COS I 



,-iEit\ 



v\) + sin( 



-_,-iE2t\ 



COS 



-^^^* + sin2^e-^^^* 






cos^sin^(e-^-^2*-e^^-^i*l 



^m) 



= Aee(t)|z/e)+^eM(i)kM) • 



(72) 



Using the above, one can compute immediately the probabilities that at time t 
the state z/e(t) is either a z/g or a v^ weak interaction eigenstate: 



P{V, -* Z/e; t) = |Aee(t)|2 = 1 " - siu^ 2d\\ - COs(^2 " ^l)t] 



P(z/e ^ Z/^; t) = |Ae^(t)| 



■sin2 2^[l-cos(E2-^i)^] 



(73) 
(74) 



Since the masses of neutrinos are small compared to the momentum, one can write 



Ei ~ \v\ + 



2 
VTi- 



Av\ 



tc^ L , 



(75) 



where L is the distance travelled by the neutrinos in a time t. Using the above, 
one can write, for instance. 



P{ue -^ z/^; L) = -sin^2e 



Am'^ 

1 — cos — r-rL 

2\p\ . 



sin 2^ sin' 



Am^L 



(76) 



where Am^ = rrig — itl1. Numerically, it turns out that 

Am'^L ^ Am2(eV^)L(m) 



4\p\ 



|p|(MeV) 



(77) 



Recapitulating, for the case of 2 neutrino species, one gets the following for- 
mula quantifying the probability that a weak interaction eigenstate neutrino (z/g) 
has oscillated to other weak interaction eigenstate neutrino (z/^) after traversing a 
distance L: 



Pii^e^ v^.]L) = sin2 2^sin2 



1.27Am2(eV^)L(m) 



|p|(MeV) 

The probability that no such oscillation took place, of course, is just 
P{ve -^ i^e;L) = 1- P{ve -^ i^f,; L) . 



(7J 



(79) 



These probabilities depend on two factors: (i) a mixing angle factor sin^ 29 and (ii) 
a kinematical factor which depends on the distance travelled, on the momentum of 
the neutrinos, as well as on the difference in the squared mass of the two neutrinos. 
Obviously, for oscillations to be important the mixing factor sin^ 29 should be of 
0(1). However, large mixing is not enough. It is also important that the kinemat- 
ical factor Am^(eV'^)L{m)/\p\(MeV) > 0(1), so that the second oscillatory factor 
in Eq. (78) can be significant. 

It is useful to develop the formalism a bit more for future use. The probability 
amplitudes Aee{t) and Ae^(t) of Eq. (72) can be recognized as matrix elements of 
a 2 X 2 matrix e"*-'^* defined by 



iHt 



where 



U 



U e-^-^diagi f/t ^ 



cos^ sin^ 
- sin 9 cos 9 



iO) 



is the mixing matrix of the weak interaction eigenstates and 

-'^diag ~ 

One has 



El 

E2 



Aee{t) 



[e-^^*]n 



Ae^,{t) 



[e-^^*]i2 



f82l 



^3) 



Using the fact that Ei = \p\ + mf/2\p\, it proves convenient to separate -f^diag into 
two different pieces: 



H. 



diag 



H + 



mf + 1712 ' 



A\p\ 



I 
1 






-1 
1 



U) 



The first piece above, because it is proportional to the unit matrix, gives an overall 
phase factor which is irrelevant for calculating the neutrino oscillation probabilities. 
Hence, effectively, one can replace 



Hrli 



diae 



H^ 






■0-3 



^5) 



In view of (85), it is convenient to define the 2x2 Hamiltonian matrix iJ^ac by 
n -U H n^-^( -cos2^ sin2^ \ 

Am 

= -^{sin2^ai - cos 2^(73} . (86) 

Then, effectively, 

AUt) = [e-^^-*]n ; A^(t) = [e-^^-*]i2 • (87) 

Just as the above coefficients describe the time evolution of a state that started 
at t = as a z/g weak interaction eigenstate [cf. Eq. (72)], one can define coefficients 

A^eit) = [e~'''-%i ; A^^it) = [e'^^-%2 (88) 

which will detail the time evolution of a state which started out at t = as a z/^ 
weak interaction eigenstate: 

|z/^(t)) = A^e(t)|z/e) + A^^{t)\iy^) . (89) 

It is easy to deduce from these considerations that the 2x2 Hamiltonian -?/vac is just 
the Hamiltonian which enters in the Schrodinger equation for |t'e(^)) and |z/^(t)): 



d 



Weit)) 



H^ 



W{t)) 
W{t)) 



(90) 



V NEUTRINO OSCILLATIONS IN MATTER 

When neutrinos propagate in matter, a subtle but important effect takes place 
which alters the ways in which neutrinos oscillate into one another. The origin 
of this effect, which is known as the MSW effect for the initials of the physicists 
who first discussed it, [II] is connected to the fact that the electron neutrinos can 
interact in matter also through charged current interactions. While all neutrino 
species have the same interactions in matter due to the neutral currents, the z/g 
weak interaction eigenstates, because of their charged current interactions, as they 
propagate in matter experience a slightly different index of refraction than the 
v^ weak interaction eigenstates (and the v^ weak interaction eigenstates). This 
different index of refraction for z/g alters the time evolution of the system from 
what it was in vacuum. 

Let us again consider the two-neutrino case. The relative index of refraction 
between z/g and z/^ is the result of the difference between the forward scattering 
amplitudes for z/g and z/^, caused by the charged current interactions of the z/g. In 
detail, one has [If] 



1 — n 



rcl 



27lNe 



Mo)L.-A{o)\ 



(91) 



where N^ is the electron density and A{0) is the forward scattering amphtude. 
The contribution of the neutral current interactions cancels in Eq. (91), while the 



iCC 



charged current contribution to A(0)|j^ ^ gives [11] 

V2GfN, 



nrel 



\P\ 



(92) 



One can use Eq. (92) and the formalism we developed at the end of the last 
section to study the evolution of neutrinos in matter. The relevant Hamiltonian 
now is 



H = i/vac + 
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(93) 



Again, because relative phases are irrelevant, we can subtract from the above a 
term proportional to the identity. This yields the following effective Hamiltonian 
describing the propagation of neutrinos in matter: 
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Because of the term proportional to the electron density in Eq. (94), in matter it 
is no longer true that the eigenstates of -f^mattcr are z/i and z/2. Calling these matter 
eigenstates u^ and z/|^, one has that 
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and 



TT _ TT rrdiag jA 

-n matter — "-^ M-" matter "-^M 



It is easy to check that -f^m'atter is given by 
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(97) 



with the mixing angle in matter 6m determined by the equation 



sin 29m = ^ ^ . (98) 

The presence of the term proportional to the electron density gives rise to inter- 
esting resonance phenomena. [12] There is a critical density A^^'''*, given by 



for which the matter mixing angle 9m becomes maximal (sin2^M -^ 1)) irrespec- 
tive of wha 
reduces to 



tive of what the vacuum mixing angle 9 is. If one is in such a medium, then -f^mafter 



rrdiag 
-'^ matter 



^ ATcrit 



Am 

sin2^a3 . (100) 
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The probability that a z/g transmutes into a v^ after traversing a distance L in 
this medium is given by Eq. (78), with two differences. First, since we are in 
a medium sin 2^ — >■ sin2^M- However, because the density is assumed to be the 
critical density, sin2^A^ -^ 1. Second, since -ffj^tter in Eq. (100) differs from Ho by 
the replacement of Am^ -^ Am'^ sin 29, such a replacement also will enter in the 
kinematical factor in the probability formula. Hence, it follows that 

(Am^ \ 

——-sin 26*1. . (101) 

4 rD / 

This formula shows that one can get full conversion of a z/g weak interaction 
eigenstate into a u^ weak interaction eigenstate, provided that the length L and 
momentum \p\ satisfy the relation 

^sin2^L = — ; n=l,2,... . (102) 

There is a second interesting limit to consider. [12] This is when the electron 
density A^e is so large that it overwhelms the other terms in -ffmatl^cr- If Gp Ne ^ 
Am^ /2^/2\p\, then one has, approximately, 

rrdiag ^F -' 'e /-. ^.qn 

-n matter — "CTS /=— • UUOj 

In this limit, it is easy to check that sin2^M -^ 0; cos29m -^ —1, so that 9m — * f- 
In this case, there are no oscillations in matter because sm29M vanishes 

-Pmatter («^e ^ l^^l^, L)\ Am2 ^ . (104) 

^ 2V2\p\Gp 



This actually is immediate also since, in this limit, i/matter itself is diagonal 



rdiag jj] Gp Ne 



-f^mattcr — U mH^^^^-^^^U j^ — (T3 -1= — . (105) 



Hence the Schrodinger equation (90), with H^^ -^ -f^matter, is diagonal and there 
can be no transitions. For future use I note that in this limit, since Om = 7i"/2, the 
Ue weak interaction eigenstate in matter coincides with the state u^: 

|z/e) = cosOmWi') + sineMW^) eM=^/2 W^) ■ (106) 



VI EVIDENCE FOR NEUTRINO MASSES 

Most experiments searching for direct evidence for neutrino masses have, up to 
now, only set limits on these masses and the associated mixing angles. However, 
there are now both strong hints, and some real evidence, that neutrino masses 
really exist coming from neutrino oscillation experiments. 

Most oscillation data is presented as an allowed region, or limits at some con- 
fidence level, in a Am^ — sin^2^ plot. That is, experimentalists find it conve- 
nient to quantify their results using the formalism discussed in the last section, 
involving oscillations among two neutrino species u^ and up The oscillation prob- 
ability formulas for Plu^ -^ i^is', L) [cf. Eq. (78)] involves both the a — (3 mix- 
ing angle Oap, as well as a kinematical factor depending on the mass squared 
difference Am^ between the mass eigenstates in the two neutrino system. Be- 
cause neutrino beams have a rather large energy spread, for Am^ ^ \p\/L the 
kinematical oscillating factor in Eq. (78) averages to 1/2. This implies that, 
in general, the sensitivity to a signal for neutrino oscillations goes down to 

Sm'^2{eaf3)mm ^ 2P(z/q ^ UfS] L) . 

A Direct Mass Measurements 

The classical way to try to infer a non-vanishing value for neutrino masses is by 
measuring /3-decay spectra near their endpoint. The presence of neutrino masses 
alters the dependence of the measured intensity /(Te) on the electron kinetic energy 
Te as it approaches the maximum energy release Q. One has [13] 

/(Te) = {Q- Te) E We^iQ " T.f - m^J ^^ . (107) 

i 

If TTiy^ = then the intensity spectrum is quadratically dependent on the energy 
release (Q — Te). If neutrinos have mass, one has a spectrum distorsion and y/l is no 
longer linear in {Q — Te), but vanishes at some value of Tg less than the maximum 
energy released Q. These distorsions are best detected in /9-decay spectra which 
have low Q values; an ideal candidate being Tritium where Q = 18.6 KeV. 



TABLE 1. Neutrino Mass Limits from ^He 
/3-decay, from Ref. [14]. 



Experiment 


"<"(eV^) 


"m,;'(eV) 


Tokyo 


-65 ±85 ±65 


< 13.1 


Los Alamos 


-147 ±68 ±41 


<9.3 


Ziirich 


-24 ±48 ±61 


< 11.7 


Livermore 


-130 ±20 ±15 


< 7.0 


Mainz 


-22 ±17 ±14 


< 5.6 


Troitsk 


1.5 ±5.9 ±3.6 


<3.9 



Tritium /3-decay experiments are sensitive to neutrino masses in the "few eV" 
range. Remarkably, most of the high precision experiments performed with tritium 
actually see an excess of events near the end-point, setting poorer limits than their 
theoretical sensitivity. [14] However, very recently, the Troitsk experiment [15] has 
been able to determine a very stringent result for the largest eigenvalue "mj^^" 
principally contributing to Tritium /5-decay: 



"™2 n 



mtj' = (1.5 ± 5.9 ± 3.6) eV ; "m^," < 3.9 eV (90% C.L.) (108) 

Table 1 gives a compilation of the existing /3-decay results in Tritium and the 
corresponding limit for "rn^y^". 

Similar, but less accurate, kinematical bounds are also known for the largest 
eigenvalues principally contributing to decays involving z/^ and z/,- weak eigenstates. 
Denoting these eigenvalues, respectively, as "m,^ " and "m,y^", one finds the follow- 
ing results. From studying vr^ -^ ^^^^^l decay at PSI [16] one has 

"m^^" = (-0.016 ± 0.028) MeV^; "m^^" < 170 KeV (90%C.L.) . (109) 

From studying the decay r — ^ z/^ + 57r at LEP [17] one has 

"m^;' < 18.2 MeV (95% C.L.) . (110) 

A different, and in some ways more interesting, limit on the neutrinos associated 
with the z/g weak interaction eigenstate comes from neutrinoless double /9-decay. 
This process, if it exists, violates lepton number. Thus it is only possible if neutrinos 
have a Majorana mass. Ordinary double /3-decay Z -^ (Z + 2) +2e~ + 2z/e conserves 
lepton number. However, in a double /9-decay processes where no neutrinos are 
emitted Z ^ {Z + 2) + 2e~, lepton number is violated. As shown schematically in 
Fig. 1, these processes can only occur if there is a neutrino-antineutrino transition 
engendered by the presence of a Majorana mass term. 

There is now a variety of measurements of ordinary double /3-decay, [18] but up 
to now there are only limits on neutrinoless double /5-decay. [19] The half-life for 
these later processes is a measure of the neutrino Majorana mass associated with 
these decays: 
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e- \ e- 

FIGURE 1. Schematic diagram which gives rise to neutrinoless double /3-decay if neutrinos have 
a Majorana mass. 
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Here the mass {mi.Jgg is given by 
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:il2) 



with rriiy. being the neutrino mass eigenstates. Note that {my^)ee vanishes if the neu- 
trinos are Dirac particles. That is, if neutrinos have only lepton number conserving 
Dirac masses. 

This point can be appreciated readily by considering, for simplicity, the case of 
one neutrino species. In this case, if these neutrinos only have a Dirac mass m/5, 
the corresponding 2x2 neutrino mass matrix M has the form 



M 



This matrix is diagonalized to 
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mr, 
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^diag = I n^^ J„ ) (114) 

by the orthogonal matrix 

U =( ^_ 2 ) • (115) 

Using Eq. (115), for this case one easily checks that {mi^Jee vanishes. One has 



TABLE 2. Bounds on ncutrinolcss double /3-decays half lives 
and associated bounds on (m^^)ee, from Ref. [14]. 



Decay 






2 


Years) 


1 


(TOi 


'Jee (eV) 


76Ge- 


-^^^se 


> 


5.7 


X 1025 


(90% 


C.L.) 


<0.2 


(90% 


C.L.) 


i28Te - 


-> i28xe 


> 


7.7 


X 10^4 


(68% 


C.L.) 


< 1.1 


(68% 


C.L.) 


i3"Te - 


^ i30Xe 


> 


5.6 


X 10^2 


(90% 


C.Li) 


<3.0 


(90% 


C.L.) 


i36xe - 


-. i36Ba 


> 


4.4 


xlO^y 


(90% 


C.L.) 


<2.3 


(90% 


C.L.) 



(m.Jee = (-^] i-mn) + (-^) (rriD) = . (116) 

Table 2 reproduces a recent compilation of experimental results on neutrinoless 
double /9-decay. [14] The most sensitive of these experiments involves the double 
/9-decay of '^^Ge to ^^Se, with a half-life limit of over 10^^ years. [20] The resulting 
bound on {m^Jee quoted is 

{rriujee < 0.2 eV (90% C.L.) . (117) 

This bound, however, has probably a factor of two uncertainty due to uncertainties 
associated with calculating the nuclear matrix elements involved in the decay. [21] 

B Cosmological Constraints 

There are some indirect constraints on neutrino masses provided by cosmology. 
The most relevant is the constraint which follows from demanding that the energy 
density in neutrinos should not overdose the Universe. Neutrinos are thermal 
relics; they decoupled from the Universe's expansion when their interaction rate 
r fell below the Universe's expansion rate H. [22] In the usual Robertson- Walker 
expanding Universe, the rate of expansion H ~ T^/Mp, where Mp ~ 10^^ GeV is 
the Planck mass and T is the Universe's temperature. Since 

T = n^{av) r^ GlT^ , (118) 

with Gf the Fermi constant, Gp ~ 10^^ GeV^^, decoupling occurs at a temperature 
Td determined by setting F ~ iJ. This gives 



1 \^''^ 



■D 



GlMp 



1 MeV . (119) 



Two cases are of interest. If neutrinos have a mass much less than To {m,j <^ To) 
they are hot relics. That is, they are relativistic at the time of decoupling. For 
hot relics, the density of neutrinos is comparable to that of photons at decoupling: 



riy r^ n^\j, . Cold relics, on the other hand, are neutrinos whose mass is much 
greater than Tn {rriy ^ Tq). In this case, at the time of decoupling the neutrino 
density n^, because of the Boltzmann factor, is much below that of the photons. 
Thus for cold relics, riij <^ n^|^ . 

In either case, one can compute the neutrino contribution to the energy density 
of the Universe. [22] This is simplest for the case of hot relics, since their number 
density essentially tracks the photon number density.^ Thus the number density of 
neutrinos now is fixed by the measured temperature of the microwave background 
radiation: 



nu = -^T^ ■ (120) 

Hence the contribution of neutrinos to the present energy density of the Universe 



IS 



Pu = Uy^m^^ . (121) 

i 

It has become conventional to normalize all densities in terms of the Universe's 
closure density pc- 






Pc = -—^ ~ 1.9 X lO^^^/i^^ ~ 1.1 X lO^/i^ — - . (122) 

SttGat cm'^ cm'^ 



In the above Hq is the Hubble constant and /i is a measure of its uncertainty. One 
finds 

Km 

Ho = lOOh 123 

sec iVlpsec 



with [24] 



Defining 



then, for hot relics, one has 



/i = 0.65 ±0.1. (124) 



fi. = - , (125) 

Pc 



92 eV h^ ■ ^ ^ 

^' Because of photon reheating at the time of recombination (and a smah statistical difference 
because neutrinos are fermions and photons are bosons), the neutrino temperature Ti, is not quite 

1 /o 

the same as the photon temperature T^. One finds T,^ — (r^) T^. [23] 
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FIGURE 2. Plot of f7j//i^ as a function of the neutrino (sum) mass. 



It follows from the above that if the sum of neutrino masses X^i "f^Vi — 30 eV, then 
neutrinos would close the Universe. Because we know that the Universe is not very 
far from closure density, if neutrinos are hot relics the sum of their masses cannot 
be much above 30 eV. Thus, although direct bounds allows for a presence of a "z/^" 
neutrino with mass "m,^ " less than 170 KeV, cosmology forbids neutrinos to have 
masses as large as that. ^ 

When neutrino masses are above Tp ~ MeV, then the simple formula given in 
Eq. (126) no longer applies. Nevertheless, it is still possible to compute Q^ taking 
into account now of the appropriate Boltzmann factor. Figure 2, adapted from [22], 
plots Qijh'^ as a function of the neutrino (sum) mass. This quantity rises linearly 
with mass up to rrij, ~ 1 MeV and then drops rather rapidly. Cosmology allows 
neutrino masses for which Qi^ < 1. So, as mentioned above, ^^my^" and ""mj,^" must 
really be well below their kinematical bounds. On the other hand, we note that our 
bounds for ^^m^J^ and {mi,^)ee lie in a cosmologically allowed region. In principle, 
cosmology also allows neutrinos to exist with masses greater than a few GeV, since 
these cold relics give Q^, < 1. However, as we discussed earlier, neutrino counting 
at LEP excludes additional neutrinos besides Ug, v^ and v^, with mass rrij^ < Mz/2. 
This exclusion region is also indicated in Fig. 2. 



^' These cosmological bounds can be avoided if the massive neutrinos were unstable and had a 
sufficiently short lifetime. [25] 
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FIGURE 3. Bounds on v^ -^ v^- oscillations for large Am^, from [27]. 



C Accelerator limits and hints for neutrino masses 



Two experiments at CERN, using z/^ beams of average neutrino energies {Ey^ ~ 
15 GeV and typical decay lengths L ~ 1 Km, put strong limits on z/^ — > v^. neutrino 
oscillations for Am^ > 10 eV^. In view of the discussion in the last subsection, 
this is a very interesting mass range to explore, since neutrinos with masses in this 
range could be of cosmo logical interest. 

These two experiments use quite different techniques to detect z/^'s. 
CHORUS [26] uses an emulsion target to try to detect the r track produced in v^ 
charged current interactions. NOMAD, [27] on the other hand, uses drift chambers 
and kinematical techniques to detect a v^ signal. The result of both experiments 
are shown in Fig. 3, along with limits obtained by some earlier experiments. The 
CHORUS and NOMAD results for Am^ > 10 eV^ exclude oscillations with mixing 
angles sin^ 2^^^ > 10~^ at the 90% C.L., improving previous limits by about a 
factor of 5. 

Because NOMAD has a very good electron identification, this experiment is also 
able to set a strong limit for z/^ -^ z/e oscillations. For ISrv? > 10 eV^, one excludes 
oscillations with mixing angles sin^ 2^^e > 2 x 10"'^ at 90% C.L. Finally, because 
the Pfj, beam at CERN has about a 1% z/g admixture, both experiments are also 
able to exclude z/g -^ z/^ oscillations for Am^ > 10 eV^, but now only for mixing 
angles sin=^ 2^er > lO-\ at 90% C.L. 

The situation regarding neutrino oscillations is much less clear cut in the region 



Am^ < 10 eV^. Here there are limits from past accelerator and reactor searches for 
oscillations and recent bounds from the KARMEN experiment. [28]. However, there 
is also some evidence for z/^ -^ v^ oscillations coming from the LSND experiment. 
[29] Let me begin by briefly detailing the results of this last experiment first. LSND 
studies neutrinos originating from pions produced at rest at the LAMPF beam stop 
by an 800 MeV proton beam. These neutrinos, which are produced in the chain 
7r+ -^ fJ'^'^ii -^ G'^'^e^fi'^fj., have average momenta in the 30 to 50 MeV range. What 
LSND looks for is the oscillation of the z/^ produced in /i"^ decay into a Ug, using 
a delayed coincidence in a target 30 meters from the beam dump. If z/^ -^ z/g 
oscillations take place, the z/g inverse /5-decay in the target (z/ep -^ e^n) produces a 
prompt photon from e~^e~ annihilation, while the produced neutron gives a delayed 
photon, as a result of the process np -^ d'-f. 

The LSND experiment observes an excess of e+7 coincidence events which, if 
interpreted as z/^ -^ z/g oscillations, give a substantial allowed region in the Am^ — 
sin^ 29 plane. However, as I mentioned above, other experiments performed in the 
past, [30] as well as the recent KARMEN experiment, [28] exclude almost all of this 
allowed region. Furthermore, new data from the KARMEN2 detector [31] which 
became available in summer 1998 appeared to exclude even the small remaining 
allowed region for LSND! 

This rather confusing situation is displayed in Fig. 4. It was discussed in some 
detail in the summary talk of Janet Conrad at the 1998 Vancouver International 
Conference on High Energy Physics. [32] As one can see from Fig. 4, a combination 
of the BNL 776 data and the Bugey reactor data only leaves the region between 
0.2 eV^ < Am^ < 4 eV^ as an "allowed" region for the LSND signal. However, 
this region is essentially excluded by the KARMEN2 data, if one uses the 90% 
C.L. bound from this experiment. However, this result itself is somewhat anoma- 
lous, since the 90% C.L. sensitivity for KARMEN2 is actually below the LSND 
signal.^ Furthermore, the LSND experiment has also looked for z/^ -^ v^ oscil- 
lations by studying z/g quasielastic scattering events and the collaboration, again, 
find an excess of events. citeLSND2 If interpreted as resulting from oscillations, 
this additional signal gives a Aw? — sin^ 26 allowed region which is consistent with 
that obtained by the z/^ — * z/g analysis. 

It is difficult to make strong statements at this stage. The best that one can say 
is that there are hints of z/^ -^ z/g oscillations in the region 0.2 eV^ < Arr? < 4 eV^, 
with rather small mixing angles sin^ 2^^e ~ 10~^. 



''"> The 90% C.L. for KARMEN2 of Fig. 4 uses data only from the initial part of their run-where 
no background events were seen, even though 3 events were expected. Additional data from 
KARMEN2 now appears to have the number of background events expected. [33] As a result, it 
looks like the full KARMEN2 results will probably be closer to the 90% C.L. sensitivity line in 
Fig. 4. 
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FIGURE 4. Summary of the experimental situation for i/^ 
Vancouver Conference, from Ref. [28]. 
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D Atmospheric Neutrino Oscillations. 



Large underground detectors, originally conceived to search for proton decay, are 
sensitive to the flux of neutrinos produced in the atmosphere. These neutrinos are 
mostly produced through the decay of pions, with the vr — >■ /z — > e chain producing 
two v^ neutrinos and antineutrinos for each z/g neutrino and antineutrino. One has 
known since the early 1990's that the observed flux of z/^'s appeared to be much 
smaller than expected, with the ratio [35] 



R 



m 



observed 






0.6 



(127) 



expected 



Although the anomalous ratio R could be the result of neutrino oscillations, 
strong evidence for neutrino oscillations only emerged in summer 1998 from the 
Super Kamiokande experiment. The SuperKamiokande collaboration [36] reported 
a pronounced zenith angle dependence for the flux of multi-GeV z/^ neutrinos, but 
no such dependence from z/g neutrinos. For neutrino energies in the multi-GeV 
range, the neutrino fluxes are not affected by geomagnetic effects in an asymmetric 
fashion. Thus one expects the observed neutrino signal to be up-down symmet- 
ric. As can be seen in Fig. 5, the SuperKamiokande data for multi-GeV z/^'s is 
clearly up-down asymmetric. There are 139 up-going u^ compared to 256 down 
going events. The observed asymmetry 
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FIGURE 5. SupcrKamiokandc results on multi-GcV i^^ events. 



7- r _ T-) s Multi-GeV 

= -0.296 ±0.048 ±0.010 (12^ 



.U + D 
is a 6(7 effect. The corresponding asymmetry for multi-GeV u, 

-jj _ r)\ Multi-GcV 



U + D 



0.036 ±0.067 ±0.020 (129) 



is quite consistent with zero. 

The SuperKamiokande collaboration [36] interprets these results as evidence for 
^11 -^ ^x oscillations, with vx some other neutrino species. This is most dramati- 
cally demonstrated in Fig. 6 where the ratio of data to Monte Carlo is plotted as a 
function of L/ Ey for both z/g and z/^ events. No LjEy dependence is seen in the z/g 
data, but the v^ data drops down to a value of 1/2 for LjEy > 10^ Km/GeV. Re- 
calling the simple 2- neutrino formula for the probability of oscillations [Eq. (78)], 
Fig. 6 suggest immediately that Am^ ~ 10~^ eV^ and that the mixing angle 9 is 
near maximal. This is confirmed by a more detailed analysis, which for z/^ -^ ux 
oscillations gives sin^ 29 = I and Am^ = 2.2 x 10"^ eV^ as the best fit point. 

It is unlikely, however, that the SuperKamiokande results are due to z/^ -^ z/g 
oscillations (that is, that ux = i^e)- First, the region in the Arri^ — sin^ 2^ plane 
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FIGURE 6. Plot of neutrino signals in the SuperKamiokande experiment as a function of L/E^ 



favored by the SuperKamiokande results, is almost totally excluded already by the 
null results of the CHOOZ reactor experiment [37] which looks at z/g oscillations into 
another neutrino species. Furthermore, the up-down ratio (129) for the z/g flux is 
more than 3a away from what one would expect if one were dealing with u^ — * Ue 
oscillations, where one expects 



U-D 
U + D 



theory 



0.205 



(130) 



If there are only three neutrino species, then most likely what is being seen in 
SuperKamiokande are u^ — > u^- oscillations. However, at this stage, it is not possible 
to rule out the possibility that ux rnay be a sterile neutrino z/^.* 

The SuperKamiokande results [36] are consistent with previous Kamiokande 
results, [38] which also had indicated a (less pronounced) zenith angle depen- 
dence of the z/^ flux. Although the Am"^ — sin^ 29 regions for SuperKamiokande 
and Kamiokande do not appear to overlap much, the 90% C.L. region of Su- 
perKamiokande is "more-signiflcant" , since the Kamiokande best flt has sin^ 26 = 
1.35. In fact, recent results presented by SuperKamiokande at DPF 99, [39] with 
more data collected, help span the gap, indicating even more clearly the consistency 
of all data with each other. 

In addition, data from other underground experiments (Soudan [40] and MACRO 
[41]) as well as other phenomena — like the flux of upward going muons [42] pro- 
duced by u^ interactions in the earth — when interpreted in a neutrino oscillation 

®^ Sterile neutrinos are, by definition, SU{2) x U{1) singlets. Because they do not couple to the 
Z, they are not excluded by the neutrino counting results from LEP. 
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FIGURE 7. Evidence for atmospheric neutrino oscillations from all experiments, from Rcf. [24]. 



framework, are totally consistent with the SuperKamiokande z/^ zenith angle re- 
sults. Fig. 7 summarizes all this information in one graph. This figure provides 
quite strong evidence in favor of neutrino masses and is probably the strongest 
evidence we have to date for physics beyond the Standard Model. 



E Solar Neutrinos 

The study of the solar neutrino flux was started in the early 1970's by Ray Davis 
and his group. [43] At present there are five different experiments which give infor- 
mation on solar neutrinos (Homestake, [44] Gallex, [45] SAGE, [46] Kamiokande 
[47] and SuperKamiokande [48]) and all five have some bearing on the issue of neu- 
trino oscillations. In fact, roughly speaking, all five experiments see approximately 
half of the expected rate, as shown in Fig. 8. However, these experiments are sen- 
sitive to different parts of the solar neutrino spectrum, because the reactions they 
use to detect solar neutrinos in their detectors have different thresholds. SAGE 
and Gallex study the reaction z/g + ^'^Ga — >■ ^^Ge + e~, which has a threshold of 
0.23 MeV. Homestake looks for the excitation of chlorine (z/g + ^^Cl — >■ ^^Ar -|- e~) 
which has a 0.8 MeV threshold. The water Cerenkov detectors, Kamiokande and 
SuperKamiokande, study elastic Uge scattering and their threshold is in the neigh- 
borhood of 6.5 MeV.^ 

It has long been felt that the observed discrepancy between the neutrino signals 
detected and the expectations of the, so called. Standard Solar Model(SSM) [49] 
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SuperKamiokande is making strong efforts to move this threshold down to 5.5 MeV. 
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FIGURE 8. Rates seen by the (liferent solar neutrino experiments, compared to the expectations 
of the Standard Solar Model. [491 



is not due to defects in this model but to the presence of some new physical phe- 
nomena. One of the principal arguments in favor of this latter solution to the solar 
neutrino puzzle, has to do with the details of the signal expected in each experi- 
ment. Because of the quite different threshold involved, each of the solar neutrino 
experiments, in fact, feels different pieces of the neutrino producing reactions in 
the solar cycle. For example, the Gallium experiments are the only ones which are 
sensitive to neutrinos originating in the pp cycle (the main solar cycle), with these 
neutrinos contributing about 50% of the expected rate. The Homestate detector 
mostly measures neutrinos from ^B, although it is also sensitive to ^Be neutrinos. 
Finally, because of their high threshold, the big water Cerenkov detectors only see 
Boron neutrinos. These circumstances make it difficult to argue for an astrophys- 
ical solution to the solar neutrino deficit. Much more natural is to imagine that 
this deficit arises as a result of neutrino oscillations. 

There are two distinct neutrino oscillation solutions to the solar neutrino prob- 
lem. Because roughly all experiments are reduced by about a factor of two from 
expectations, it is possible to fit the data by using vacuum neutrino oscillations 
^'e — ^ ^x- Clearly, for this fit one must appeal to large mixing angles and assume 
a tiny Am^. Since E^ ~ MeV and the earth-sun distance L ~ lO^"*^ m, typically 
Am^ ~ 10^^^ eV^. However, because the Homestake result is only about 30% of 
the predicted value, one has to fine-tune the parameters, so that only a few "just 
so" regions are favored. [50] A recent "just so" fit by Bahcall, Krastev, and Smirnov 
[51] is shown in Fig. 9. 

In my opinion, much more interesting that the above "solution" is the possibility 
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FIGURE 9. "Just so" solar neutrino fit, from [51]. 



that the solar neutrino results are a reflection of matter induced oscillations (the 
MSW effect we discuss in Section V). In the sun, the electron density to a good 
approximation, can be characterized by an exponential profile function [52] 



Neir) = Ne{0)e «o 



:i31) 



The central density A^e(O) — 10^^ cm ^ ~ 10^^ (eV) is rather high and for appro- 
priate values of Am^ and sin^ 29 can exceed the critical MSW density. For instance, 
for 9 small, Am^ ~ 10"^ eV^ and \p\ ~ 3 MeV 
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(132) 



It follows from Eq. (132) that, for these parameters, z/g's produced in the core of 
the sum (where Ne{o) 3> iV™*) as they radiate outward go through a region with 
Ne ~ A^g'''* and can oscillate to z/^'s (or other neutrino types) without paying a 
mixing angle penalty, since sin 26*^ I /vcrit -^ 1- 

The actual calculation of what happens in the sun is rather complicated, [12] 
since the density N^ changes along the neutrino trajectory. Since Ng = Ne{t), the 
matter Hamiltonian of Eq. (94) is now time dependent: 
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Although one can diagonahze this Hamiltonian, the resulting mixing angles and 
energies will be time dependent: 
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Because of this time dependence, it is no longer true that the states \i^i^{t)) and 
|z/^(t)) are actual eigenstates. In fact, transitions can occur between these states. 
A simple calculation [12] shows that the states \i^^(t)) obey a coupled Schrodinger 
equation 
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then transiting between |z/{^(t)) and |^'|^(t)) will be relatively unimportant and 
one has an adiabatic situation. For an exponential density profile, Eq. (137) is 
satisfied at A^"'* provided that [12] 



Am^ sin^ 20 
2 cos 26 



> 2 X 10-^ (eV)" . 



(138) 



For the adiabatic case, one can use our discussion of matter oscillations to give 
a qualitative picture of how the MSW mechanism could work in the sun. Because 
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at the solar core A'e(o) 3> A^e"*; according to Eq. (106) |z/e) 
we are assuming adiabaticity, as the neutrinos diffuse out of the core of the sun, 
the state |t'|^(o)) will evolve into \i^^it)). That is, there are no transitions in the 
sun. Thus, when the neutrinos exit the sun, the state |^'|^(tsurfacc)) will just simply 
become |z/2). Because 



(^^6^2) = sin6' , 



(139) 



it follows that, in this case. 
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A more careful analysis shows that there are actually two MSW solutions, one 
adiabatic and one non-adiabatic, [53] both having Am"^ ~ 10~^ eV^. The adiabatic 
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solar neutrino data, from Ref. [51]. 



sin 2Q plane favored by the MSW explanations of the 



solution has large mixing angles sin^ 1Q 
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1. Hence, according to Eq. (140), 
indeed, roughly half the flux is lost. The non- 
adiabatic solution has sin^ 26* ~ 5 x 10^^. Furthermore, a rather large range in 
the IStt? — sin^ Id plane is eliminated by the absence of a day /night effect, which 
would be a sign of matter oscillations in the earth. The favored MSW regions are 
depicted in Fig. 10. 

I want to close this brief discussion of solar neutrinos, and in particular of the 
MSW explanation of the solar data, by making a more quantitative remark. To fit 
the solar neutrino data using the MSW effect requires that the probability P(z/e — * 
Ve,\ L) have considerable energy dependence. The required energy dependence is 
shown in Fig. 11. I indicate also in this figure at what energies the neutrinos 
produced in the various solar reactions are effective. One sees from Fig. 11 that 
essentially all pp neutrinos survive, the Berylium neutrinos disappear and the flux 
of Boron neutrinos is roughly halved. This reconciles nicely with what is seen in 
the data, as detailed in Table 3. 
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Table 3. Summary of neutrino flux observations and expectations of solar 

neutrino experiments 



3^01 (13% Be; 80% B) 
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^^Ga (51% pp; 15% Be; 12% B) 

(77.5 ± 7.7) SNU Ref.[45] [predicted SSM : (129 ± 8) SNU] 
(66.6 ±8.0) SNURef.[46] 
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VII THEORETICAL IMPLICATIONS 

I summarize in Fig. 12 where one is at present on the issue of neutrino masses and 
mixings. This figure collects together all the neutrino oscillation evidence, as well as 
the hints for oscillations, which we have at the moment. As is clear from the figure 
there are three regions suggested in the Am^ — sin^ 29 plane. The strongest evidence 
is that for atmospheric neutrino oscillations coming from the SuperKamiokande 
zenith angle data. Here the suggested parameters are (Am^) ~ 3 x 10~^ eV^, 
sin^26' ~ 1 with u^ -^ ux {i^x 7^ t'e)- Solar neutrinos also are strongly suggestive 
of oscillations. Interpreting the data this way leads to (Am^) ~ 10~^ eV^, with 
sin^2^ ~ 1 or sin^2^ ~ 5 x 10"^, for MSW u^ -^ ux oscillations, or Am^ ~ 
10~^^ eV^ and sin^2^ ~ 1, for "just-so" z/g ~* t^x oscillations. The weakest hint 
for oscillations probably is that of LSND, because of other contrary evidence. At 
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FIGURE 12. Summary of evidence of, and hints for, neutrino oscillations. [24] 



any rate, the suggested region here is Am? ~ 5 x 10 ^ eV^ and sin^ 20 ~ 10 ^ for 



l^f: 



u^ oscillations. 



Besides neutrino oscillation phenomena, there is really no other direct evidence 
for neutrino masses. However, both /9-decay and neutrinoless double /9-decay put 
rather strong bounds on neutrino masses connected to z/g. From /3-decay the largest 
neutrino mass obeys the bound: "mj,^" < 3.9 eV (90% C.L.). Double /3-decay, 
bounds the Majorana mass {m^Jee even more strongly: {my^)ee < 0.2 eV (90% 
C.L.). These bounds are interesting, since they are close to the kind of neutrino 
masses which could have substantial cosmological influence. Using the central value 
for the Hubble parameter [c.f. Eq.(124)], I note that Y^i'^Vi = 30 eV, 6 eV, 2 eV 
correspond, respectively, to neutrinos closing the Universe, to neutrinos being 20% 
of the dark matter in the Universe (assuming VLm — 1), and to neutrinos being 
20% of the dark matter, with VLm — 0.3. The last case is, perhaps, the one that 
is most cosmologically realistic. [55] This stresses the the importance of continuing 
the search for neutrino masses in the eV range. 

Concentrating only on the SuperKamiokande evidence for neutrino masses al- 



ready has important implications. Taking 

Am^ =ml-mlr^3x 10"^ eV^ (141) 

gives one already a lower bound on some neutrino mass: m^ > 5 x 10^^ eV. 
This mass value, in turn, gives a lower bound for the cosmological contribution of 
neutrinos 

Although this number is far from that needed for closure of the Universe, I note 
that the Qi, of Eq. (142) is comparable to the contribution of luminous matter to 
the energy density of the Universe [56] 

^luminous ~ (3 - 7) X 10"=^ . (143) 

So, from SuperKamiokande we learn that the neutrino contribution to the energy 
density of the Universe is the same as, in the words of Carl Sagan, that of "billions 
and billions of stars" ! 

For particle physics, a value of rris ~ 5 x 10~^ eV is also quite interesting. If we 
use either the simple see-saw formula of Eq. (55) or the GUT relation (66), the 
identification 



rs^ 



[("^d)3]" 

ms={ ™s 0^2™^ (144) 

rriT — ^^^^ 

give comparable values for ms and A: 

m5 ~ A ~ 10^^ GeV . (145) 

Of course, these values are only justifiable in specific models where one has a bit 
more control of other constants, which are taken above all to be of 0(1). 

If one goes beyond the SuperKamiokande data, then many theoretical scenarios 
emerge. Unfortunately, in general, these scenarios mostly reflect the prejudices one 
has regarding the data. Nevertheless, it is useful to briefly discuss two differing 
broad theory scenario. In the first scenario, one assumes that all hints for oscilla- 
tions seen are true. In the second, one disregards some oscillation hints. In most 
cases, the discarded data is that of LSND. 

If one believes all hints for neutrino oscillations, since there are three different 
Am^ involved, the neutrino mass matrix M necessarily is a 4 x 4 matrix. ^° To get 
a 4 X 4 neutrino matrix one adds to the usual three neutrinos a sterile neutrino Ug. 
Most 4 neutrino models attempt to fit all data, since this was after all the reason 
for introducing the fourth neutrino. The most promising scenario [58] has two pairs 



^"^ There have been attempts to "stretch" some of the data, so that all hints can be accounted 
for with only two different Am^. These attempts [57] seem rather forced to me. 



of quasi-Dirac neutrinos split by a small mass difference. ^^ The heaviest pair {1712 
and m^) have masses of order 0.5 eV and contribute Q,, ~ 0.03 to cosmology. The 
atmospheric neutrino oscillations involve this pair, so that Am^^ ~ 3 x 10"'^ eV^. 
The second pair {mi and 7714) are much lighter, with mass around 10~^ — 10~^ eV. 
Their mass difference Am^4 ~ 10~^ eV^ is what enters in solar neutrino oscilla- 
tions. The LSND result is explained as an oscillation between the light pair and 



the heavy pair, with Am 



~ m. 



0.6 eV . In this scheme, the solar neutrino 



oscillations involve oscillations of u^. to a sterile neutrino, while the atmospheric 
neutrino oscillation is u^ -^ v^ and LSND v^. — > z/^. Although this scheme works 
phenomenologically, theoretically it is difficult to get light sterile neutrinos almost 
degenerate with ordinary neutrinos. 

Different patterns arise if one is prepared to disregard some of the neutrino 
oscillation limits. If one disregards, in particular, the results from LSND then 
the CHOOZ bound, [37] and the quite different mass squared differences involved 
in atmospheric and solar neutrino oscillations, suggest a very simple 3-neutrino 
mixing matrix. CHOOZ suggests that ^13 ~ 0°. On the other hand, atmospheric 
neutrino oscillations suggest ^23 — 45°. Finally, depending on what solar neutrino 
oscillation solution one picks, the angle 612 can either be large or small. Thus, 
neglecting possible CP violating phases, the neutrino mixing matrix looks like [59] 
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(146) 



where C12 = cos ^12; S12 = sin ^12. However, even with U of the above form, there 
are many open questions to answer. For instance, is maximal mixing [612 — 45") 
allowed? Are nearly degenerate neutrino masses (mi ~ m2 — m^ ~ 0.5 eV) 
allowed? What neutrino mass matrix gives rise to this particular mixing matrix? 

These questions cannot really be answered in a straightforwad manner, without 
making some more assumptions. There is really a lot of freedom. Given U and some 
assumptions for the neutrino mass spectrum {mi} then one can deduce a neutrino 
mass matrix M. However, recall from our discussion of the see-saw mechanism, that 
M itself depends on both the neutrino Dirac mass mo and on the right-handed 
neutrino mass matrix ms [cf. Eq. (55)]. Thus, to make progress, even "knowing" 
M one has to make some assumptions on itid (or ms) to learn something further. 
For instance, one could use GUTs which naturally ties the matrix mD in Eq. (55) 
to the M-quark mass matrix. [60] 

Some authors have preferred to focus on some simple structure for the 3x3 
matrix M . [61] Two of these are particularly appealing. The first of these has 

^^> A quasi-Dirac neutrino pair reduces to a Dirac neutrino, as the mass difference between tfie 
pair vanislies. 



a total degeneracy for the neutrinos, the other is Dirac-hke with an additional 
massless neutrino. In the first pattern 



M = m 1 . (147) 




In this case, J2i itT'u, = 3m, so that cosmology impose a bound on m, depending on 
what one believes Q^j is. However, in the degenerate case, one has also that 

{'m,Jee = J2u!^m,^=m. (148) 

i 

The double /3-decay bound then tell us that m < 0.2 eV. If we push m to its upper 
bound, however, it is difficult to see what perturbation can then give Am 
3 X 10-3 eV^; Am^^i^, ~ 10"^ eVl 

The second simple pattern for neutrino masses has [62] 
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which has a degenerate pair and a zero eigenvalue. Note that this pattern conserves 
Lf, — L^ — Lt-. Since for this mass matrix Aml^^^^ = m^, it follows that m ~ 
5 X 10^^ eV. So in this case, neutrinos do not contribute much to the energy 
density of the Universe {Q,y ^ 3 x 10"'^). To get solar neutrino oscillations one 
has to introduce some perturbation on the mass matrix (149) that will split the 
massive degenerate states and give Am^g ~ 10"^ eV^. 

VIII FUTURE EXPERIMENTS 

It seems pretty clear that progress in understanding what is going on in the 
neutrino sector can only come from further data. Fortunately, new data will be 
forthcoming in all the relevant Arn^ regions. I want to end these lectures by briefly 
discussing these future experiments. 

A Solar Neutrinos 

SuperKamiokande will continue to take data in years to come, thus refining their 
present measurements of solar neutrinos. Furthermore, a real effort is taking place 
to lower the neutrino energy threshold further so as to be able to study the shape 
dependence of the signal as a function of E^,. In addition to this continuing effort, 
relatively soon two other experiments will be coming on line which have considerable 
promise. The first of these is SNO (the Sudbury Neutrino Observatory [63]) which 



uses a Kiloton of D2O. The advantage of having heavy water is that it allows SNO 
to study simultaneously both charged current and neutral current processes. The 
charged current process 

h'e + d^'e~+p + p, (150) 

like all charged current processes, is sensitive to whether oscillations have occurred 
or not. The neutral current disintegration of the deuteron, on the other hand, is 
insensitive to oscillations since it is the same for all neutrino species I'x'^'^ 

ux + d^ux+p + n. (151) 

Comparison of the rates for the two neutrino reactions(150) and (151) should help 
rule out possible astrophysical explanations for the solar neutrino puzzle. The SNO 
detector should begin taking data in 1999. 

The second solar neutrino experiment of interest is Borexino. [64] This experi- 
ment is presently under construction at the Gran Sasso Laboratory and should be 
ready for data taking in 2001. Borexino uses 300 tons of scintillator, which has a 
relatively low threshold (-Ethr > 340 KeV). As a result, Borexino should be par- 
ticularly sensitive to the E^, = 862 GeV neutrino line coming from ^Be. Recalling 
Fig. 11, one sees that if the MSW explanation is correct, the solar neutrino signal 
in Borexino should be significantly below the theoretical expectations. Indeed, if 
there are no solar oscillations, Borexino is supposed to detect about 50 events/day, 
while if the MSW explanation is true, this number should go down to about 10 
events/day. 

B Atmospheric Neutrinos 

Here again SuperKamiokande will continue to integrate data with time. However, 
the Am^ region, will also be probed more directly by using neutrino beams from 
accelerators. Three such long baseline experiments are in different stages of readi- 
ness. K2K, which uses a neutrino beam from KEK, aimed at SuperKamiokande 
250 Km away, should shortly be operational. [65] MINOS, [66] in the Soudan Mine, 
is under construction and will be the target of a dedicated neutrino beam from 
Fermilab, 730 Km away. First data should become available around 2001-2002. 
Finally, a variety of proposals exist for experiments in the Gran Sasso Laboratory, 
which is 740 Km from CERN, to become targets of neutrino beams from CERN. 

The main advantage that these long baseline experiments have over Su- 
perKamiokande is that the neutrino beam used is well characterized, both in energy 
and in its time structure. Furthermore, these beams also have a higher intensity. 
So many possible systematic effects will be under better control. In the case of 
the higher energy Fermilab and CERN beams, it may also be possible to directly 
search and detect z/t-'s, if these neutrinos are produced in the oscillations. 



^^•^ This is only true for neutrinos whose neutral couplings to the Z are universal. It does not 
apply to sterile neutrinos. 



C LSND Region 

The Am^ — sin^ 29 region identified by the LSND experiment as potentially in- 
teresting also will be explored further. At Fermilab, there is an approved experi- 
ment, Mini BooNE, [67] which will run around 2001-2002, which should be about 
a factor of five more sensitive than LSND in a comparable kinematical region. 
With this sensitivity, it should be quite clear whether v^ -^ z/g oscillations with 
Am^ ~ (0.1 — 1) eV^ exist or not. 
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APPENDIX A: DIRAC AND MAJORANA MASSES 

To understand how Dirac and Majorana masses can arise, it is useful to review 
here some of the properties of the spinor representations of the Lorentz group. 
The Lorentz group, besides the well known vector and tensor representations has 
also spinor representations. It turns out that there are two inequivalent spinor 
representations. It is out of these two-dimensional spinors that one builds up the 
usual four-dimensional Dirac spinor ip. 

Under a Lorentz transformation, a vector field V^ has the well known transfor- 
mation 

V^' -^ V'^ = A^^V , (1) 

where the 4-dimensional representation matrices A obey the pseudo-orthogonality 
conditions 

%u = KVapA^. , (2) 

involving the metric tensor 



%,y 



(3) 



Besides vector representations, the Lorentz group has two inequivalent spinor 
representation. The corresponding 2-dimensional Weyl spinors are conventionally 



denoted by C,a and ^a, known as undotted and dotted spinors, respectively. Under 
Lorentz transformations they transform as 

Ca^C = MX (4) 

ia^C = M:% . (5) 

The 2x2 matrices M and M*, with det M = det M* = 1, provide inequivalent 
representation of SL{2, C). Obviously, from the above it follows that ^ ~ ^*. 

One can establish a relationship between the 2x2 matrices M and the 4x4 
matrices A, since the vector field V^ transforms as V ~ ^ ® ^. For these purposes, 
it is useful to define a set of four matrices a^ = (1, a), with a being the usual Pauli 
matrices. The 2x2 matrix 

V = a^^v.uV = a>^V^ (6) 

under a Lorentz transformation transforms as 

V ^V' = MVM^ = af'V^ . (7) 

Using Eq. (Al), it follows that 

<a; = My^,,Mf . (8) 

Because det M = 1, the analogue of the scalar product for vectors V^rifj^iyV^ = 
V^V^, for the spinors C, and C, leads to the following Lorentz scalars: 

UX = U' ; UX = W (9) 

where e"* = — e*" and e^^ = 1. Similarly, just as the contraction of a covariant and 
contravariant metric tensor gives the identity [rj^prjP'^ = 6'^], one can define 2x2 
antisymmetric e-matrices, eat, which obey 

eace^' = ^l ■ (10) 

It follows that ei2 = —1. 

The usual 4-component Dirac spinor ip is made up of a dotted and an undotted 
Weyl spinor: 

* - (I" ) (11) 

In this, so called, Weyl-basis the Dirac 7-matrices 7^^, which obey the anticommu- 
tation relations {7^, 7^^} = —21]'^'^, take the form 



Here cr^ = (1, —a), so that in this basis 



7« ^ ? i h .- 1, ^: (13) 



^10/' ' \ -a' 

and 



75 = ^lVlV = " 1 • (14) 



It follows from the above that C,a and x"" are chiral projections of ip: 

^L = ^(1 - 75)^ = ( ) ' ^ = V^^(l + ^5) = (0 C) (15) 

^R = ^(1 + 75)^ =(;['«)' ^ = V^^(l - 75) = ir* 0) . (16) 

Using these equations, it is easy to see that the Dirac mass term connects ^ with 
X- Specifically, one has 

/:Dirac = -mol^^R + ^^l) = -mD(Cx" + ^^L) ■ (17) 

Recall, however, that dotted spinors are related to the complex conjugate of an 
undotted spinor. Choosing a phase convention where 

C = L; x: = Xa, (18) 

one can write the Dirac mass term simply as 

CBir.c = -mniiar + x''Q ■ (19) 

In view of Eq. (A9) this term is obviously Lorentz invariant. However, Lorentz in- 
variance does not require one to have two distinct Weyl spinors C, and x to construct 
a mass term. Majorana masses, basically, make use of this "simpler" option. 

One can define a 4-component Majorana spinor in terms of the Weyl spinor C, 
and its complex conjugate C,'- 

^M = ( ^: ) . (20) 

Because ^" = ^"*, effectively ipM has only one independent helicity projection. One 
can choose this projection to be, say, {4'm)l'- 

(V'm)l = ^(1 - 75)^M = ( o" ) ; "(^^ = ^^(1 + ^5) = (0 L) ■ (21) 



One can construct {iPm)r by using the charge conjugate matrix C. In the Weyl 
basis C is given by 
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(^m)r = (r 0) = ie%, 0) = (ae^a 0) = {iPM)lC . 
That is, {iPm)r is the charge conjugate of {iPm)l (c-f- Eq. (36)): 

[(^m)l]^ = (^m)r • 



(23) 



(24) 



Because of Eq. (A24) it follows that the Majorana spinor ipM obeys a constraint. 
It is self-conjugate: 



^a 

r 



V'M = I r 1 = (^m)l + (^m)r = (^m)l 

Hence, 

The Majorana mass term 

-^Majorana = — -'m M'4' M'4' M 

involves a product of ^ with itself and ^ with itself 
1 



[(V^m)l]^ . (25) 



(26) 



^Majorana = --^M ( (V^m)l(V^m)r + (V^m)r(V^m) 



(27) 



:^M(ear + rea). (28) 



Eq. (A28) can also be written purely in terms of {iPm)^ by using the charge 
conjugation matrix C. Using Eq. (A23) one has also 



-^Majorana = -""^M ((V^m)lC(V^m)l + (V'm)lC'(V'm)l 



(29) 



Equally well one can write this mass term entirely as a function of {iPm)r- One 
finds 



1 



-^Majorana = -^ ((V^m)rC(V^m)r + (V^m)rC(V'm)r 



T 



(30) 
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